Spin correlations in the algebraic spin liquid - implications for high Tc superconductors 
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We propose that underdoped high Tc superconductors are described by an algebraic spin liquid 
(ASL) at high energies, which undergoes a spin-charge recombination transition at low energies. 
The spin correlation in the ASL is calculated via its effective theory - a system of massless Dirac 
fermions coupled to a U{1) gauge field. We find that without fine tuning any parameters the gauge 
interaction strongly enhances the staggered spin correlation even in the presence of a large single 
particle pseudo-gap. This allows us to show that the ASL plus spin-charge recombination picture 
can explain many highly unusual properties of underdoped high Tc superconductors. 
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I. INTRODUCTION 



One of the intriguing questions about the cuprate su- 
perconductors - to many possibly the key to understand- 
ing superconductivity - is the role played by spin corre- 
lations in these materials. By now it is well established 
that the insulating parent compound of the copper oxide 
superconductors is well described as a two dimensional 
Heisenberg Antiferromagnet in the temperature regime 
above the three dimensional Neel ordering temperature 
(Tat ~ SOOX). Between Tn and a temperature of about 
600 K the materials seem to follow the renormalized. clas- 
sical results of Chakravarty, Halperin and Nelsontl. At 
higher temperatures the system appears to cross over 
to quajutum scaling behavior with a dynamical exponent 
z — iB. On doping with holes, away from stoichiom- 
etry, these insulating compounds develop into high Tc 
superconductors even for very low hole concentration of 
order 5%. For such low doping the low energy Hilbert 
space should still be dominated by spin fluctuations - as 
a memory effect of the undoped parent compound. The 
question that needs to be addressed is the peculiar in- 
terplay of short range antiferromagnetic correlations as a 
remnant of the ordered Neel state at zero doping compet- 
ing with spin singlet formation present in the supercon- 
ducting state. Depending on the importance attached to 
those two spin correlations we can distinguish two schools 
of thought that have strongly influenced the high Tc re- 
search. 

The first takes the strong antiferromagnetic (AF) 
spin fluctuations seriously as providing the "mechanism" 
for superconductivity. Assuming the normal state is a 
Fermi liquid, the antiferromagnetic spin fluctuation^ were 
shown to lead to a d-wave superconducting stateE The 
spin fluctuations can be rnpdeled phcnomenologically a la 
Millis, Monien and FinesQ in the context of a nearly an- 
tiferromagnetic Fermi liquid (NAFL). The main degrees 
of freedom are spin and charge carrying quasiparticles 
coupled to their own spin fluctuations. The strength of 
the spin fluctuations is modeled by a susceptibility which 
is strongly peaked near the commensurate antiferromag- 
netic wavevector Q = {■k / a^ir / a). This strong enhance- 



ment is achieved via an RPA in the form 



x(q,t^) = 



xo(q,'^) 



1- J(q,r)Xo(q,t^) 



as proposed by Monthoux and Pines (xo(q, tij) is calcu- 
lated by taking spin fluctuations mediated interactions 
betweoa the quasiparticles into account). In a series of 
papersH Monthoux and Pines have given what they call 
a "proof of concept" for spin-fluctuation driven pairing 
for the cuprates near optimal doping. They study the 
competition of spin-density wave (SDW) order with su- 
perconducting pairing by solving two-dimensional Eliash- 
berg equations for spin-fluctuation induced interactions 
between the quasiparticles. These approaches could ac- 
count for both the d-wave pairing symmetry and the 
high Tc's seen in the cuprates. One interesting question 
is of course the interplay of SDW and superconductiv- 
ity as instabilities of the NAFL. Monthoux and Pines 
concluded that superconductivity always prevails as the 
dominant instability. With the realization that the quasi- 
particles are not too close to a SDW instability they how- 
ever needed a strong fine tuning of their Stoner enhance- 
ment factor J(q) as a function of both momentum depen- 
dence (strongly peaked at Q = (7r/a,7r/a)) and strength 
(>^(Q)xo(Q, 0) ~ 0.97) to account in a "self-consistent" 
manner for their empirical input spin spectrum. 

The discussion so far was centered on the optimally 
doped case, in this paper we shall however be concerned 
with the underdoped cuprates where the spin behavior 
shows the peculiar competition of antiferromagnetic or- 
der and singlet formation in a particularly striking way 
as is evidenced by various spin-pseudo-gaps seen in NMR 
and Neutron scattering. Let us briefly mention the ap- 
proaches to this regime which have emphasized the anti- 
ferromagnetic order and the concomitant spin wave exci- 
tations as the dominant degrees of freedom. In an exten- 
sive study of the 2c? antifercpmagnetic Heisenberg model 
Chubukov, Sachdev and Ya3 have identified the spin cor- 
relations in the quantum critical regime at the transition 
from Neel order to paramagnetism. In their paper they 
also argue for the appearance of quantum critical scaling 
away from zero doping by incorporating the effect of the 
doped charge degrees of freedom into a reduction of the 
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spin stiffness (finite in the ordered Neel state of the 2d 
Heisenberg model) Sokol and PineaZl have similarly ar- 
gued for quantum critical behavior and its cross over to 
a quantum disordered regime driven by the frustration 
of antiferromagnetism via the hole motion as the physics 
behind the strange spin correlations in the pseudo-gap 
regime. They postulate a susceptibility which is dom- 
inated by the effect of a finite coherence length in the 
quantum disordered regime and the appearance of the 
corresponding full energy gap for spin excitations. 

As we shall argue below this phenomenology qualita- 
tively agrees with what we obtain within the slave bo- 
son approach to the < J-model, where however in con- 
trast to the above, spin singlet formation (or more pre- 
cisely, the formation of the d- wave/staggered- flux ghase 
- a state with a particular quantum order UlCnOlnBu) is 
the driving force behind all of the strange phenomenology 
in the underdoped samples. The slave boson approach is 
a microscopic theory for high Tc superconductors which 
takes the on-site repulsion U and the resulting charge gap 
(~ 2eV) as the largest scale in the theory. It naturally re- 
produces the Mott insulator phase at half filling and low 
charge carrier density at low doping. Furthermore this 
framework allows us to calculate various line-shapes when 
the AF-correlation is short ranged. The line-shapes can 
be quite different from those obtained from phenomeno- 
logical models. 

The slave boson approach is strongly tied to the strong 
coupling phenomenology incorporated in the Hubbard or 
i J-model and was guided by Anderson'sEij exciting pro- 
posal of a spin liquid as a realization of the strongly cor- 
related Mott insulator in the parent compound of the 
high Tc superconductors. More specifically, Anderson 
proposed that the cuprate physics could be understood 
in terms of doping holes into a state which consists of 
preformed spin singlet dimers on nearest neighbor bonds 
whose quantum fluctuations lead to uniform strength 
AF-correlations on all nearest neighbor bonds (termed 
resonating valence bond (RVB) state). The RVB state 
also has short range antiferromagnetic spin correlations 
(no true long range order) due to the singlet formation 
and hence comprises a spin liquid phase. However, due 
to the success of the two dimensional Heisenberg antifer- 
romagnet in describing the insulating parent compound 
of the cuprates this idea of a spin liquid ground state was 
partly rejected. 

Nevertheless the RVB picture has a lot of appeal. Let 
us mention briefly the rational behind the opinion that 
the spin singlet formation dominates antiferromagnetic 
correlations in the underdoped cuprates. The slave boson 
approach treats the underdoped cuprates as doped Mott 
insulators - the key experimental fact of high Tc super- 
conductors. It also naturally incorporates the spin sin- 
glet formation and predicted, ahead of experiments, the 
pseudo-gap metallic phasellJ in the underdoped cuprates 
and the superconducting d-wave order (preferred over 
s-wave by strong-. Coulomb correlations present in the 
Mott insulator) .E3 The pseudo-gap metallic phase is a 



new state of matter and it is very rare in the history of 
condensed matter physics that a new state of matter was 
predicted before its experimental observation. An impor- 
tant effect of the preformed spin singlets, present in the 
RVB picture, is the fact that the spin on the doped holes 
can become an excitation on top of the RVB collective 
state whereas the charge remains tied to the empty site. 
This led to the notion of separate spin and charge car- 
rying excitations called spinous and holons respectively. 
Superconductivity then appears as a consequence of es- 
tablishing coherence among the doped holes which natu- 
rally incorporates a low (proportional to x the hole dop- 
ing concentration) superfluid density present in the un- 
derdoped cuprates. This picture of the spin fluctuations 
is quite different from the NAFL where (1) a metallic 
state with no pseudo-gap is assumed (which does not 
apply to the underdoped regime) and (2) strong antifer- 
romagnetic spin fluctuations couple to spin and charge 
carrying quasiparticles. It is the exchange of those anti- 
ferromagnetic fluctuations which gives rise to the pairing 
necessary for superconductivity. 

However, underdoped cuprates have a very puzzling 
property which seems hard to explain using the spin liq- 
uid approach (or other approaches). As the doping is low- 
ered, both the pseudo-gap and the AF correlation in the 
normal state increase. Naively, one expects the pseudo- 
gap and the AF correlations to work against each other. 
That is the larger the pseudo-gap, the stronger the spin- 
singlet formation and the weaker the AF correlations. 
One result coming out of the research reported here is 
a connection of the quantum scaling regime of the un- 
doped compounds at high temperature with Anderson's 
RVB picture via the algebraic spin liquid (ASL) phase. 
We find that the ASL can explain both the pseudo-gap 
and the strong AF correlations in the underdoped high 
Tc superconductors. Due to the U{1) gauge fluctuations 
in the ASL phase, the AF spin fiuctuations in the ASL 
are as strong as those of a nested Fermi surface, de- 
spite the pseudo-gap. The doped holes suppress the U{1) 
gauge fluctuations which explains the decrease of the AF 
spin fluctuations with increasing doping even though the 
pseudo-gap is smaller at higher doping. We would like 
to stress that the ASL approach to the AF fluctuations 
is distinct from RPA-like approaches. Due to the small 
density of states in the pseudo-gap phase, one needs to 
fine tune a large coupling constant within the RPA to fit 
the experimental data for underdoped samples. Such a 
fine tuning is not necessary in the ASL approach. Thus 
the U{1) gauge structure (and the associated spin-charge 
separation) in the ASL lie at the heart of the strange phe- 
nomenology of the_underdoped cuprates in the so called 
pseudo-gap phaseE^I. 

The starting point of the slave boson approach is a 
microscopic lattice model. Among the popular models, 
the tJ-GtJ ( Generalized-i J) model seems particularly 
promising as a description for the low doping regime 
where the competition between delocalization energy t 
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and the spin fluctuations J becomes manifest. 

The form of the tJ Hamihonian can be justified starting 
from the Hubbard model in the limit of strong on-site 
Coulomb repulsion energy which in the cuprates leads to 
an insulating charge transfer gap of 2eV . As the Coulomb 
energy is the largest energy scale in the problem it is nat- 
ural to treat the kinetic energy as a perturbation which 
to lowest order in leads to the tJ Hamiltonian with 
the single site Hilbert space restricted to the three states 
spin up/down and empty. For zero doping the tJ Hamil- 
tonian reduces to the Heisenberg model with antiferro- 
magnetic exchange coupling (virtual hopping on top of 
strong onsite Coulomb repulsion naturally leads to anti- 
ferromagnetic exchange). 

The organization of the paper is as follows; In the 
next two sections (|l|lll) we review the ideas underlying 
the ASL and discuss the spin correlations in this phase 
sketching the calculations involved (details can be found 
in the appendix). Following this we give a brief discus- 
sion on-how the recently introduced concept of quantum 
orderSEl affects the current results (IV). Section (|^ ad- 
dresses the problem of what happens to the ASL physics 
at low energies where the interactions between spinous 
and holons become strong. This will be followed up by 
a discussion of possible implications for the cuprates and 
a short comparison with other approaches to the same 
region of the phase diagram (Vl). The final section (VII) 
will summarize the main results and discuss open prob- 
lems for further study. 



II. SU{2) SLAVE BOSON APPROACH AND 
ALGEBRAIC SPIN LIQUID 

Taking the charge gap as the largest energy scale, we 
have argued that the natural starting point to the de- 
scription of underdoped cuprates is a strongly correlated 
model. The question of spin correlations in the cuprates 
ties into the question of how to reconcile local moment 
magnetism with itinerant electron spin density fluctua- 
tions which is a longstanding problem in condensed mat- 
ter physics. There seems to be consensus in the com- 
munity on how best to describe the two extreme lim- 
its of the cuprate phase diagram - the Mott insulator 
at zero doping as a Heisenberg antiferromagnet and the 
heavily doped regime in terms of itinerant electron mag- 
netism. Of course as always the most interesting regime 
is the one lying in between the two limits which shows 
the antiferromagnetic-singlet dichotomy described above. 
For the lightly doped regime one could argue that clearly 
a perturbation about the Neel ordered state with frustra- 
tion due to holes, a la Sokol and PinesEl, makes the most 
sense. However, this approach does not specify the na- 
ture of the quantum spin liquid in the disordered phase. 



Understanding the properties of the quantum spin liquid 
is very important since it is the quantum spin liquid that 
controls the pseudo-gap metallic state in the underdoped 
cuprates. Another drawback with the NAFL approach 
to the spin and charge fluctuations is the fact that it 
cannot explain naturally the small charge carrier density 
without breaking translation symmetry. 

The picture of spin fluctuations followed here is the one 
of preformed singlets (RVB) which incorporates the spin 
correlations in the spin disordered phase (such as the su- 
perconducting state and the pseudo-gap normal state) in 
a natural way by paying the prize of throwing away a lot 
of antiferromagnetic correlations. The aim of the current 
paper is to show a way of reconciling the singlet for- 
mation with the strong dynamic antiferromagnetic cor- 
relations seen in NMR and Inelastic Neutron Scattering 
(INS) with the help of a gauge theory based on a slave 
boson decoupling of the < J-model. 

As we discussed briefly above, the (generalized) tJ- 
model can be shown to be an excellent description of the 
low energy physics embodied in the Hubbard model un- 
der the condition that the single site Hilbert space is con- 
straint such that double occupation is forbidden. On an 
operative level this problem can be attacked by introduc- 
ing slave particles that change the constraint cj^Co-(i) < 1 
(where c„{i) is the physical electron destruction operator 
of spin a on site i) into 'Saflfai'i) + b'^b{i) — 1 where 
the b{i) is a bosonic operator - the holon - that keeps 
track of the empty sites and carries the charge of the 
physical hole, whereas the spin of the hole is carried by 
fa - the fermionic spinon. This is achieved by writing 
Cff = fab^ and should be read as an equality in the con- 
straint Hilbert space. With the introduction of the slave 
bosons we have however introduced a redundancy in the 
description of our physical problem which is a U{1) phase 
rotation of / and b that leaves the physical electron op- 
erator invariant. As the system evolves under the tJ 
Hamiltonian (|l|) this phase will strongly fluctuate as a 
function of space and time. It is this phase degree of 
freedom that corresponds to U{1) gauge fluctuations and 
makes the description of the constraint problem of spin + 
hopping in terms of fermionic and bosonic operators pos- 
sible. The above program can be summarized under the 
name of U{1) slave boson theory and was mplemented 
widely in the early |tla,ys of high Tc researchllil. Early on 
it was also realizedEj that the Heisenberg model in the 
fermionic representation has an extra SU(2) invariance 
organizing the spinous into doublets 



hi 



fli 



(2) 



In the usual U(\) slave boson approach this invari- 
ance was lost on introducing holes. More receally Wen 
and Lee introduced a slave boson formulationEjilj which 
maintains this SU{2) structure away from half filling. 
This was achieved via the introduction of slave boson 
doublets and representing the physical electron operator 
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more accurately called UlCnOln phase. Quantum order 
J and its PSG characterization is a very important concept 

C|i = --j=h\il>^i — --j=(b\^f^i + fo2i/|i) (3) for our discussion (See section IV). 

ex. = ^hU^, = 7l(^i^^i^ - ^2./!.) gpj^ SUSCEPTIBILITY 



with 



the doublet of bosonic fields keeping track of the doped 
holes. 

The slave boson approach to the tJ Hamiltonian 
then is ta. u.se-.this decoupling and perform a mean-field 
analysisllja'Ej. This has led to a mean-field phase dia- 
gram as a function of hole doping x and temperature T 
which is in qualitative agreement with the phase diagram 
of the cuprates. Within the mean-field description how- 
ever the gauge freedom is treated only on the average - 
i.e. replaced by a static configuration of phases which 
in turn determine the band-structure of the spinous and 
holons. However as is usually the case with the identi- 
fication of phases via a mean-field decoupling we need 
to consider fluctuations about the mean-fields to deter- 
mine their stability. In particular the gauge field - un- 
constrained by any dynamics - will affect the physical 
properties in each phase drastically. 

In the present paper we are concerned with the pseudo- 
gap regime of the cuprates which - within the U{\) for- 
mulation - was identified as the d-wave paired state for 
the spinons. Within the SU{2) approach this phase can 
also be described as the so called staggered-Flux (sF) 
phase, without explicit fermion pairing. As it turns out 
(seetZl for more details) the SU (2) formulation allows for 
a more straightforward identification of low-lying (mass- 
less) gauge modes and it was shown that the sF phase 
breaks the SU{2) gauge structure down to U{1). This 
massless U{1) mode was missed in the early discussions 
(within the U{1) slave boson approach) on the pseudo- 
gap phase and plays a crucial role in the pseudo-gap 
phase of the underdoped cuprates. It is responsible for 
the emergence of the ASL. In the next section we shall 
discuss the physical spin correlations in the ASL. 

Before closing this section, we would like to remark 
that the ASL (the sF phase) is only one of many possible 
symmetric spin liquids. All symmetric spin liquids have 
the same symmetry, and hencewe cannot distinguish dif- 
ferent symmetric spin liquids by their symnrietries. The 
concept of quantum order was introducedB'D to distin- 
guish the different internal structures present in the sym- 
metric spin liquids. For the symmetric spin liquids con- 
structed within the SU{2) slave boson approach, one can 
use the projective symmetry group (PSG) to character- 
ize their different quantum orders. One finds that the sF 
phase is described by the particular PSG with the name 
UlCnOln and it is one of an infinite number of possible 
symmetric C/(l) spin liquids. Thus the sF phase can be 



Our starting point is the sF (or UlCnOln) state in 
the SU (2) mean-field phase diagram where the effective 
degrees of freedom are spinons and holons coupled to 
a massless U{1) gauge field. In order to analyze this 
problem we have mapped the lattice effective theory for 
the sF state (at zero doping) onto a continuum tiieory of 
massless Dirac spinors coupled to a gauge field Jl3 whose 
Euclidean action reads 



S 



N 



(4) 



where Va-.o — 1 and TV = 2, but in the following we will 
treat N as an arbitrary integer. In general Va,i ^ 
However, for simplicity we will assume Va^i — 1 here. The 
Fermi field is a 4 x 1 spinor which describes lattice 
spinons with momenta near (±7r/2, ±7r/2). The 4x4 
7^ matrices form a representation of the Dirac algebra 
{7^1, Iv] — 2(5^1/ {[I, = 0, 1, 2) and are taken to be 
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-CT2 



(5) 
(6) 



with CT^ the Pauli matrices. Finally note that = 
^'3^70- The dynamics for the J7(l) gauge field arises solely 
due to the screening by bosons and fermions, both of 
which carry gauge charge. In the low doping limit, how- 
ever, jwe will only include the screening by the fermion 
fields which yields 



exp 



d3 
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(27r)3"^^^^^^^ 



H„ 



(7) 



By simple power counting we can see that the above po- 
larization makes the gauge field a marginal perturbation 
at the free spinon fixed point. Importantly however we 
should note that since the conserved current (that cou- 
ples to a^) cannot have any anomalous dimension, this 
interaction is an exact marginal perturbation protected 
by current conservation. 



A. Uniform correlations 

Let us now discuss how the gauge fluctuations affect 
spin correlations near momentum transfer q = (0,0). 
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FIG. 1: Non-zero leading corrections to the staggered 
spin correlation function. The x denotes the vertex which is 
the 4x4 unit matrix in the case of interest. 

The expression for the uniform spin correlation reads (see 
Appendix ^) 

(^^(^)^„"(0)> - / j^e^'-'(S:{q)S-i-q)) (8) 

{S+{q)S-i-q)) = -\j ^Tr[7oG(p07oG(p~g)] 

where (• ■ ■) denotes the expectation value with respect to 
theory (||) . From this expression we see that the uniform 
spin correlation is proportional to IIoo, the polarization 
operator of the spinons. Hence it cannot be strongly af- 
fected by the massless gauge field as we argued above via 
current conservation. This was shown in an explicit cal- 
culation byr-Chen, Fisher and Wu in the context of the 
FQH effectc^L Contrasting to this as was suggested by 
Kim and Leelij we would expect the gauge fluctuations to 
strongly affect the staggered spin correlations which are 
not protected by current |Gonservation and restore anti- 
ferromagnetic correlation^^! which have been largely lost 
in the mean-field singlet state as we shall see next. 

B. Staggered correlations 

In the following we will present fluctuation corrections 
to the antiferromagnetic spin-spin correlations at order 



1/N . The expression for the staggered correlation is ob- 
tained in the form (see Appendix |A|) 

{St{q)S-{-q)) = -\j -^T,[tG{p)tG{p-m 

At the mean-field level this is simply 

{St{q)S-{-q))^ = -\ j ^Tr[lGo(p)lGo(p-9)] 

where Go{p) — and the vertices are the 4x4 unit 
matrices denoted 1. Decoupling the denominator in the 
usual way via Feynman parameters we obtain 

{Stiq)S-iq))o^-^^^^ (10) 

Note that q is measured with respect to (tt, tt) in units 
of the lattice spacing. From expression ([To|), we see that 
the antiferromagnetic correlations have been largely lost. 
The reason for this can be traced back to the spin singlet 
formation in the mean-field sF (UlCnOln) phase. This 
motivates the inclusion of gauge fluctuations. At order 
1 /N we have the three non- vanishing diagrams depicted 
in Fig. 01 . In order to calculate the contribution of these 
diagrams we note that unlike the single spinon spectral 
function the density density correlation is gauge invariant 
and hence we can choose to work in the Landau gauge 
where the gauge propagator reads 

DM<1) = -^{S,.-^-^) (11) 

and we have the following expressions for the diagrams 
shown 
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^ i ijf^ i ^ i jo,'' i £) (ij) 

i{p + k),Y i{p + k + q)sY i{k + q)plP ikaj" 



[1(G)] =- 



d'^k 



df^q 



{2ttY J {2ttY 



Tr 



(12) 



(13) 



After performing the trace and integrating over fc in c? = 3, where the integrals are convergent (see Appendix]^), 
we arrive at 



2p- q + f 



p-q 



\P\ 



[1(G)] 



\q\^\p + q\ \p\q^\p + q\ \q\\ 
1 r d'^q I 2 \p\ f + 2p-q 4|pl q^f + 2p^ 

N J (2^ [f^W~ JWW+^ ~ f\p + q\ ' p-q<f\p\\p + q] 

I 



(14) 
(15) 



In order to proceed we need to regularize the above in- 



tegrals. Because of the dot-product appearing in the de- 
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nominator of the last term in (|15|) it is hard to use dimen- 
sional regularization. Hence we have set d — S and intro- 
duced an upper momentum cutoff A. Thus performing 
the final integrals over q we can extract the ^oij-divergent 
term in the form 

§iA)] + 1(5)] + 0(0] = _^|p-lZn(^) 
I 



After combining this with the mean-field result ( |lO| ) 
and continuing to real frequencies (see Appendix ^) we 
obtain 



lmx(t^,q) = X (^^,q) = Im(5'+(w,q)S'^ (-t^,-q)> 



32 



(16) 



where C's is a constant depending on the physics at the 
lattice scale. In the limmit N 00 this reduces to the 
mean-field result ( |lO| ) 

From (|l^) it is clear that the gauge fluctuations have 
reduced the mean- field exponent. If we boldly set = 2 
which is the physically relevant case we find > 1/2 sig- 
naling an antiferromagnetic instability. This is consistent 
with an earlier numerical result obtained in Ref. In 
Fig.|^ we plot the imaginary part of the spin suscepti- 
bility at q = = Qaf = (tt, tt) 

This result is quite natural in the light of what has 
been said so far. The gauge fluctuations arise from the 
constraint of no double occupancy. The dynamics for the 
gauge field - in the very low doping limit assumed here - is 
solely due to virtual spinons. It is the low energy spinon- 
antispinon pairs having a nesting condition for scatter- 
ing between the nodes (separated by Qaf) m the sF 
(UlCnOln) phase mean-field spectrum which are respon- 
sible for the antiferromagnetic enhancement in the ASL 
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FIG. 2: Imaginary part of the spin susceptibility at Qaf- 
Note the divergence at small cu 



phase. Let us emphasize again that this enhancement 
of spin correlations at Qaf = (tt, tt) over the uniform 
component is protected, just hke the gapless U{1) gauge 
fluctuations and the gapless spinons, by the UlCnOln 
quantum order 

We have thus established algebraic behavior for the 
staggered spin correlations. It is such an algebraic cor- 
relation that leads to the name ASL that we assigned to 
the phase. 



C. Correlations near (tt, 0) 

For this momentum transfer the spinons get scattered 
between the two different types of nodes in the sF spec- 
trum which yields the following expression for the corre- 
sponding correlations 

(S(i.o,(?ii=r.,o)(-<B) = (") 

Having worked hard to obtain the anomalous dimen- 
sion of the staggered spin operator by brute force we will 
here resort to standard field theory renormalization ap- 
paratus which allows for a more economical derivation 
of the anomalous dimensions of composite operators. In 
order to obtain the anomalous dimension of the spin op- 
erator which is a fermion bilinear we need to obtain its 
wavefunction renormalization in the form 

Zs = Z^Zt (18) 

where F is the relevant vertex and as we have seen 
above depends on the momentum transfer (F — 

7o,l, 'q^ ^ near (0, 0), (tt, tt) and (tt, 0) respec- 

tively) . 
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FIG. 3: Fermion-Fermion vertex 



The spinon wavefunction renormalization is obtained 
from the self energy in the usual way and evaluates in 
the Landau gauge to 



= 1 



4 1 T{^) 



(19) 



7V37r2 (Ap){3-d)/2 

where M is the renormalization scale and we have used 
dimensional regularization. 

To obtain the vertex renormalization we need to eval- 
uate the divergent part of the 1-loop diagram depicted in 
Fig.Q which reads 



{2n) 



where F is the relevant vertex and is given by (11). 
Since the vertex correction diverges only logarithmically 
we can simply evaluate the divergent part by setting all 
the incoming momenta equal to 0. Hence the above sim- 
plifies to 



N J {2TrY 



;5 



where we now concentrate on the (tt, 0) vertex. 
This can be easily evaluated with the help of 

7er(^,0) = (1 ~ 2^e2)r(^,o)7e 



(2^13' 



1 



47r2 3 - d 



and results in 



4 1 r(^) 



(20) 



which together with equations ( p^ ) and (|l^) gives 

^5,„,„, = 1 (21) 

Hence the spin operator at (tt, 0) doesn't pick up any 
anomalous dimension from the gauge field interactions. 
It is not hard to check within the above outlined calcu- 
lation that for the other two vertices corresponding to 
the uniform and staggered spin operator, we recover the 
results discussed in the previous sections which gives a 
nice check on the brute force calculation. 



IV. PERTURBATIVE STABILITY OF THE sF 
(UlCnOln) PHASE 

An important technical detail was hidden in the above 
calculation. In this section, we will expose this issue. 

The sF phase contains 2 families of 4-component mass- 
less Dirac fermions, which couple to massless U{1) gauge 
fluctuations. At the perturbative level, the interaction 
can generate many possible counter terms (in this case 
the self-energy terms) which can 

(A) generate a mass for the Dirac fermions, 

(B) generate a finite chemical potential for the Dirac 
fermions and change the Fermi points into Fermi pock- 
ets, 

(C) shift the crystal momenta of the Dirac fermions from 
(±7r/2, ±7r/2) to some other values. 

In our perturbative calculation in the last section, we 
have assumed that non of the above counter terms are 
generated by the interaction. 

In a generic interacting theory, we used to believe that 
all the counter terms that are consistent with symme- 
tries will be generated. All of the above three types of 
counter terms are consistent with the underlying lattice 
symmetries and are hence allowed by the symmetries. In 
fact they do appear in our calculation of the continuum 
theory as cut-off dependent terms. It appears that the 
results in the last section are incorrect since the impor- 
tant counter terms are ignored. Those counter terms, if 
present, will drastically alter our previous results. 

As was stressed in this paper and in Ref.|| and Refj^, 
the sF phase, as a quantum spin liquid, is not only char- 
acterized by its symmetry but also by its quantum order. 
The quantum order in the sF phase is characterized by 
a particular PSG - UlCnOln. The mean-field ansatz of 
the sF phase is invariant under the transformations of the 
UlCnOln PSG. As pointed out in Ref.| and Ref. |, the 
perturbative fluctuations around the mean-fleld ground 
state may deform the mean-fleld ansatz. Those deforma- 
tions correspond to the counter terms. However, in our 
case, it is incorrect to use symmetries to determine which 
counter terms are allowed. One should use the PSG to 
determine the allowed counter terms. This is because 
perturbative fluctuations can only deform the ansatz in 
such a way that the deformed ansatz remains invariant 
under the same PSG. It was shownO that the above three 
counter terms are forbidden by the UlCnOln PSG. This 
is why we can drop them in our calculation. 

When we calculate the effects of interactions in a con- 
tinuum theory, we introduce a high energy/short distance 
cut-off. This cut-off destroys the structure of the under- 
lying quantum order. Thus it is not surprising that the 
counter terms forbidden by the quantum order show up in 
continuum theory as cut-off dependent terms. To restore 
the quantum order that existed in the underlying lat- 
tice theory, we can simply drop all the forbidden counter 
terms in our calculation within the continuum theory. 
We see that the PSG and the concept of quantum order 
play an important role even in calculations within con- 
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tinuum theories. It is the understanding of the PSG and 
quantum order that makes sensible calculations in the 
continuum limit possible. In fact, the theory of quan- 
tum order is partly motivated by the above issue of the 
counter terms. 

In summary, the perturbative fluctuations around the 
mean-field sF (UlCnOln) state cannot change the quan- 
tum order and cannot generate jenergy gaps for the 
U{1) gauge field and the spinonsEI. Thus the gapless 
[/(I) gauge fiuctuations and the gapless spinous are pro- 
tected by the UlCnOln quantum order present in the 
sF (UlCnOln) phase. They remain gapless even when 
the interaction is finite all the way down to zero en- 
ergy. The interacting gapless excitations lead to many 
unusual properties of underdoped cuprates, such as the 
non-Fermi liquid behavior of the ncmnal metallic state, 
the broad electron spectral functiorO, and the diverging 
AF spin fluctuations in the presence of the pseudo-gap. 
The UlCnOln quantum order in the sF phase not only 
protects the gapless excitations, it also protects the mo- 
mentum of the gapless excitations. For example the spin- 
1 gapless excitations can only appear near k = (0,0), 
(tt, tt), (0,7r), and (tt, 0). 

We would like to point out that although the ASjUcan 
be a stable quantum phase in the large N limit p'E3 for 
the real N = 2 case, non-perturbative instanton effects 
cause an instability. Thus at low energies the ASL will 
change into some other state, such as the d-wave super- 
conducting state, the AF state, stripe states, or even 
a Z2 spin liquid state. Since the UlCnOln quantum 
order in the sF phase requires the spin-1 gapless exci- 
tations to appear at k = (0, 0), (tt, tt), (tt, 0), (0, tt), the 
shift of the low frpqiiprify iTPiitrpij <ipattpriTicr peak ob- 
served in experimentsS'U'OCjCZlo'Ej indicates a transi- 
tion from the ASL to some other state at low temper- 
atures. Studying how the momentum of spin-1 gapless 
excitations shifts away from (7r,7r), (tt, 0) and (0,7r) will 
allow us to experimentally identify the low temperature 
phase. In the last section, we have studied the spin fluctu- 
ations in the ASL. In the following section, we will study 
the spin fluctuations in the low temperature phase. 



THE FATE OF THE ASL 



In section III, we showed how the ASL physics recon- 



ciles the pseudo-gap formation (ie the condensation into 
spin singlets within the RVB picture) with enhanced dy- 
namical antiferromagnetic fluctuations. At lower ener- 
gies the ASL is unstable - and an important question is 
how the ASL evolves in the underdoped regime into the 
superconducting state as we reduce temperature. 

In the following we are going to address this question 
via first analyzing the effect of the opening of a gap in 
the gauge fluctuations on the staggered spin correlations. 
Intuitively we should expect that for energies below the 
mass gap the spin correlations should be given by mean- 
field correlations - which is indeed the case. Having es- 



tablished the profound effects of this gap formation we 
take a more careful look at the change of the mean-field 
correlations on going from the sF (UlCnOln) phase into 
the fermion pairing state which will allow us to address 
the question of incommensurate spin fluctuations seen at 
low frequencies and temperatures in the cuprates. Let us 
now proceed to consider the effect of giving a mass to the 
gauge fluctuations. In the context of the spin correlations 
considered here, we would thus expect the destruction of 
the antiferromagnetic enhancement below the mass gap 
and a resurfacing of the singlet character of the correla- 
tions in the mean-field sF (UlCnOln) phase. 

To implement the mass gap formation we take 
the foUowiftg phenomenological form for the gauge 
propagatorll3 



The calculation then goes through as above (with one 
subtlety in the analytic continuation discussed in the Ap- 
pendix ^ resulting in the following expression for the 
staggered correlation 



Im(5+(c.,q)57(u;,q)) 
= 6(^2 - q2)e [m^ + q2 - C, 



(22) 



2v 



where 



32 

In Fig.Q we show the resulting spectra with m = 
20meV. Fig.||^ depicts the spectrum at the antiferro- 
magnetic ordering wavevector Qaf- As expected the 
opening of the massgap in the gauge fluctuations has re- 
stored the mean-field result and consequently suppressed 
the antiferromagnetic enhancement seen in Fig.B. 

From this analyzes we see that the spin correlations 
at low frequencies are described by the underlying mean- 
field . In the above calculation we have assumed that the 
only effect of the destruction of the massless U{1) gauge 
structure in the sF phase is the formation of the t/(l) 
mass gap which leaves the underlying mean-field correla- 
tions in tact at low energies. On a more profound level 
however the massless nature of the U{1) gauge structure 
in the sF siate is a manifestation of the UlCnOln quan- 
tum ordeiTJ present in this state (see introduction and 
section iV). Thus it is not possible to simply break the 
gauge structure without affecting the underlying mean- 
field ansatz as we have assumed above. If we reverse the 
above logic, this implies that the low frequency spin cor- 
relations which reflect the underlying mean-fleld ansatz 
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Energy [meV] 



FIG. 4: Scan of the imaginary part of the spin susceptibihty in 
arbitrary units for q = — tt/IO . . . tt/IO and lu = . . . GOmeV . 
We have chosen m = 20meV in this plot 
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FIG. 5: Imaginary part of the spin susceptibihty at Qaf 
where m = 20meV notice the recovery of the linear mean- 
field result below the massgap 



can be used as an experimental probe into the corre- 
sponding quantum order/gauge structure accompanying 
the mean-field correlations. 

The ASL described by the sF (UlCnOln) state may 
change into several different states at low temperatures. 
The low energy spin correlation can have different behav- 
iors in those low temperature states. In the following, we 
will discuss one particular low temperature state to gain 
some intuition on the behavior of the low energy spin 
correlation in the low temperature phases. 

Within the SU{2) slave boson theory, the mean- 
field phase which is energetically favored over the sF 
(UlCOln) phase at-low temperatures is the so called d- 
wave pairing statetZl. In order to keep this section reason- 
ably self contained we would like to outline some of the 
steps leading to the mean-field correlations discussed in 
the following. As was mentioned above within the SU(2) 
slave boson formulation the physical electron operator is 



represented as: 



(23) 



' Id 



where the following SU(2) doublets were introduced 



hi 



hi = 



The ip^i, ipii are the two fermion fields representing the 
destruction of a spin up and spin down on site i respec- 
tively and hi is the doublet of bosonic fields keeping track 
of the doped holes. Putting this representation into the 
t-J Hamiltonian 



H 



E 

(y) 



1 



[J{S, ■ Sj - -n,nj) - t{cl^c„j + h.c.)] 



yields on performing a Hubbard-Stratonovich transfor- 
mation to the appropriatCpiosonic bond variables the 
following partition functional 



J DhDh^ D^}^ D'il>DaoDUe ^ 



(24) 



L 



<i3> 



m{\i'l^r'^.^ + hWh, 



a 



(25) 



The oo fluctuations incorporate the projection to the 
space of SU(2) singlets within the above representatiou. 
for the electron operators (p3[). Furthermore note thatEHl 
J = 3J/8, t — t/2 and the matrix Uij in the form 



-X^ 

A*; 



Ay 

Xij 



contains the Hubbard-Stratonovich flelds which classify 
the part of the phase diagram we are looking at. The 
mean-field phase diagram is found by minimizing the free 
energy for a given number of particles with respect to the 
bond variables Uij. 

The d-wave pairing state can be represented as 

Ui^i+n = -T^Y-hr^A 



(26) 



The key difference to the sF (UlCnOln) phase is the 
appearance of a finite Og which acts as a chemical poten- 
tial for the spinous. Without this term the above ansatz 
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is S'C/(2)-gauge equivalent to the sF (UlCnOln) ansatzO. 
It is shown in RefJ^ that the appearance of such a chem- 
ical potential term is not consistent with the quantum 
order of the sF (UlCnOln) phase and hence signals the 
presence of a different quantum order accompanying the 
d-wave pairing state. It is precisely a nonzero Og which 
gives rise to an anomalous fermion fermion pairing which 
when combined with holon condensation leads to the d- 
wave superconducting state. As we shall see in the fol- 



lowing, Oq is also responsible for the incommensurate spin 
response at the lowest frequencies shifting the peak away 
from Qaf = (ti", tt). 



Taking the above ansatz equation (26) we can calcu- 
late the spin-spin correlations at mean- field level (see Ap- 
pendix |E) to obtain for the imaginary part of the spin 
correlation near Qaf = ('"',71') 



lm{S+iu, Qaf + 1^)5- (-w, -Qaf - k)> 



9 [uj^ - {vfki - 2alf - v^kfl ^ uj'^ - [vfki - 2alY - v^k^ + 
9 [uj^ - {vfk2 + 2alf - vlkl] ^uj-^ -{vfk2 + 2alY -vlkl + 
9 [uj^ - {vfk2 - 2alf - v^kj] - {vfk2 - 2a3)2 - v^kf 



ki 



kx k" 



V2 



fc2 



V2 



Vf=2V2aJx V2 = 2V2aJA 
I 



(27) 



& 

V 1 



(0= 0.1 J 



ky 




kx 



-1 -1 




FIG. 6: Wave vector scan of x" {kx, ky,u! = O.IJ). Notice the 
incommensurate pattern at[4«w energy with amplitude peaks 
shifted along the diagonalsEiJ (0, 0) corresponds to Qaf and 
kx, ky are measured in units of 1/a with a the lattice spacing 



In Fig.Q we plot the resulting spectrum for a fre- 
quency Lo < LUc = (^) ''o where the intensity is 
peaked around four points shifted diagonally away from 
Qaf- In contrast to this in Fig.Q we depict a scan for 
uj > iOc the spectrum is peaked at four points shifted hor- 
izontally (tt ± 5, tt), (tt, tt ± S) where the ridges - growing 
out perpendicular to the diagonals - overlap. S is related 
to Oq via Sa = ^ with a the lattice spacing. It should 
be noted that in the uj ^ the peaks are always located 



FIG. 7: Wave vector scan of x"ikx,ky,uj = O.IJ), kx,ky are 
measured in units of 1/a. With agchosen to give 5a — Q.l r.l.u. 

(27r = 1 r.l.u. )@Elthe value oiujc is set by lj^ = ^ (^^^ SaJ. 

To compare with experiments we have chosen the values — 



1/70 which gives lOc = 0.2J. Notice how the overlapping 
ridges lead to the horizontal incommensuration. 



along the diagonal, as expected, for in this regime the 
response is dominated by creation of particle hole pairs 
connecting the Fermi points. Interestingly for intermedi- 
ate energies the peak intensity in the spectrum is shifted 
horizontally away from Qaf even without one dimen- 
sional phenomenology. 
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VI. APPLICATION TO HIGH 

Before launching into a discussion of how the above 
mentioned physics embodied in the ASL might shed some 
hght on the cuprate spin physics we need to briefly digress 
and estabhsh some nomenclature. As we shall see in the 
following each experimental community has their own set 
of pseudo/spin-pseudo-gaps and it is easy to get lost in 
the confusing terminology. Here we adopt the following 
convention: 

pseudo-gap phase denotes the part in the cuprate phase 
diagram where angle resolved photoemission experiments 
observe a large gap in the single particle spectrum at 

The term spin- pseudo- gap will be used in conjunction 
with the spectrum at the antiferromagnetic wavevector 
Qaf = (tt, tt) to denote the frequency range of reduced 
spectral weight below the peak above Tc (which on cool- 
ing below Tc shifts to higher frequencies and becomes the 
"resonance" - see discussion below). 

Over the years NMR and inelastic neutron scattering 
(INS) have painted a very interesting picture for the spin 
correlations in the cuprates starting with the discovery 
of the by now famous AlmeV "resonance" peak by J. 
Rossat-Mignod et a&B in the superconducting state of 
near optimally doped YBCO. In addition to this reso- 
nance mode which appears exclusively in the supercon- 
ducting state there is now evidence for precursory pil- 
ing up of spectral weight in the energy region of the 
"resonance" even in the normal state of underdoped 
cuprates. The energy at which this enhancement occurs 
in the normal state is sometimes referred to as the "spin- 
pseudo-gap" energy by INS workera^j (hence our conven- 
tion). From the perspective of spin fluctuations INS is 
clearly the most powerful probe giving both dynamical 
and wavevector dependent information about the spec- 
trum. 

NMR on the other hand probes particular regions in 
wavevector space with Knight shift probing the uniform 
susceptibility - cx xi'i = 0, a; ^ 0) - and spin lattice re- 
laxation rate probing ~ i™x(qc„.o,") ^^^^ where qc„ = 
Qaf and qo = are the characteristic wavevectors de- 
termining the spin lattice relaxation at the planar Cu and 
O sites respectively. In addition to the above mentioned 
"spin-pseudo-gap" , NMR has identified two characteris- 
tic temperatures. One is the reduction of the Knight shift 
below a temperature which we will denote by T/v which 
indicates the loss of spectrum for x(q — 0,uj ^ 0) be- 
low this characteristic temperature. Furthermore there 
is a second lower scale which shows up in the spin lat- 
tice relaxation rate for Cu as a peak in (T) which we 
shaU call Tt^. Notice that Tt^ < T/v but both of these 
temperatures are above the superconducting transition 
temperature Tc in underdoped samples. From the per- 
spective of mean-field theory it is impossible to explain 
these observationsEll. Within the mean-field picture the 
reduction of the Knight shift below T/v can be related to 
the reduction of spin fluctuations on condensing into the 
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FIG. 8: Imaginary part of the spin susceptibility at Qaf 
taken from Sternlieb et al. Note the temperature dependence 
of the Imx for a fixed cu converges to a universal function of 
T/lu at large T. 



spin singlet sector identified as the sF (UlCnOln) phase 
above. However there is then no way to explain the dif- 
ference in the spin lattice relaxation rates observed for 
planar C„ and O. Whereas ^^j^{T)c^ increases with de- 
creasing temperature between T/v > T > Tt-^, jhp{T)o 
reduces monotonically below T/v. 

Let us now discuss these interesting phenomena from 
the point of view of the ASL physics described in the. 
previous part. As has been suggested by Kim and LeeE2l 
the enhancement of C„ over O relaxations is expected 
to depend on the inclusion of the gauge fluctuations. 
The rational for this suggestion arises from the fact 
that Cu — NMR probes q ~ QAP-Auctuations which 
are strongly affected by the gauge fluctuations whereas 
O — NMR derives its main contributions from q ^ 
which is protected from gauge fluctuations by current 
conservation. We have shown above that this enhance- 
ment indeed takes place with ImxQAF,i.j^o (see Fig.|^) 
being strongly increased over the mean-field behavior 
which determines the spin fluctuations at q = 0. Given 
the different T- dependence for the spin lattice relaxation 
rates for temperatures T/v > T > Tt^ it is thus natural 
to identify this regime_with the ASL physics. 

Sternlieb et aS have performed tempera- 
ture dependent neutron-scattering experiments on 
YBa2Cu30e,e{Tc — 53K) and found a universal 
dependence of 



Imx(^) = J Imx(q, w)dq cx Imx(QAF,t^) 



(28) 



on the ratio uj/T above a temperature T ~ lOOK. Their 
data (see Fig.|8|) suggests that the temperature at which 
this universal scaling appears increases with decreasing 
energy. In particular for uj = 5meV the onset temper- 
ature is T lOOA'. Below this characteristic temper- 
ature (which again is well above Tc) their data shows a 
decrease in spectrum which they associate with the open- 
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FIG. 9: Temperature dependence of l/T\T{Cu) and 
l/TiT{0) in yB02C„306.63. The dotted-dashed line shows 
the temperature dependence of the static susceptibiUty. The 
data is taken fromcl who quotes Takigawa's resuhs 
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FIG. 10: Phase diagram implied by NMR experiments. 



ing of a spin-pseudo-gap below tOg ~ \QmeV . They also 
point out that their data and the fact that the onset 
temperature of scaling increases with decreasing probe- 
frequency are consistent with the anomalous behavior of 
T^(^)c„ which probes a; ^ and increases down to 
T = Tt, ~ 150K 

As we argued above it is exactly the regime for lo > m 
at r ~ or T > Tm at that we associate with the 
ASL physics. Under the assumption of perfect scaling we 
can thus account for the difference in (T) between C„ 
and O quite naturally. From 



— «E 



Imx(q,cj) 



(29) 



where for C„: (^^ '-^ Qaf) due to the form factor 
whereas O probes the uniform spin fluctuations. In the 
ASL we can for C„ approximate X]q ^^^^^TT^ U^o ^ 
imx(QAF,t>j) 1^^^ which via scaling results in 
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(30) 



whereas i^^{T)o falls monotonically as T is reduced 
following the uniform susceptibility (see Fig.|^). As 
2i/ — 1 ^ this is consistent with the weak tempera- 
ture dependence seen by Sternlieb et al Fig.||] at large 
T/w. 

The ASL gets destroyed below the massgap energy 
scale which we associate with a temperature Tm- Be- 
low this temperature scale both the C„ and the O spin 
lattice relaxation rates fall monotonically following the 
sF-singlet correlations as T is reduced. This is seen in 
experiments (see Fig.[||) 

Thus we associate the temperature Tt^ ^ I'aQK where 
7Yt(-^)C" shows the peak with Tm the temperature be- 
low which the massless U{1) gauge structure (and thus 



the underlying quantum order) becomes unstable and un- 
dergoes a "transition" (see discussion at the end of this 
section) 

To summarize the above comparison between the ASL 
and NMR experiments, we present the phase diagram 
Fig.jl^. Below the pseudo-gap scale T^g, which we asso- 
ciate with T/v, the oxygen 1/TiT starts to decrease due 
to the opening of the pseudo-gap associated with the spin 
singlet formation. However, the copper l/TiT keeps in- 
creasing as temperature decreases, despite the small sin- 
gle particle density of states in the pseudo-gap regime. 
This strange behavior can be explained very well by the 
ASL due to the diverging AF spin fluctuations at low en- 
ergies even in the presence of the pseudo-gap. Below the 
temperature Tm, the U{1) gauge field starts to gain a gap 
and the enhancement in the AF spin fluctuations ceases 
to exist. This causes the copper 1/TiT to decrease with 
decreasing temperature following the expected mean- field 
behavior. Thus the ASL described by the sF state (or the 
UlCnOln state) appears between Tm (associated with the 
experimental temperature scale Tt-^) and Tpg (associated 
with the experimental scale T/v). Below Tm, the ASL will 
change into another state whose nature will be discussed 
later. 

Let us now discuss the INS results which in their own 
right paint a very interesting picture of the spin fluctu- 
ations in the cuprates. First we consider the spectra at 
the antiferromagnetic ordering wavevector as a function 
of temperature and energy. Already well above the su- 
perconducting transition in the normal state, INS data 
show a marked increase in spectral weight at a finite char- 
acteristic frequency E^pg and a reduced spectral weight 
for uj < Espg ( the spin- pseudo-gap see Fig.[0 bottom 
curve and compare Fig.@). O n cooling below Tc this 
gap increases in size and the corresponding peak shifts 
to higher energies as depicted in Fig. [ll|] where we have 
added the two dashed lines to the data taken by Fong et 
aB^ to guide the readers eye. From the perspective of the 
ASL liquid the "normal state" above Tc corresponds to 
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Peak Susceptibility at (0.5,0.5) 




FIG. 11: Energy dependence of peak intensitjr^t Qaf as 
a function of temperature taken from Fong et a£3. Note the 
shift in the peak position on decreasing the temperature below 
Tc — 67K (the dashed hues were added by the current authors 
as a guide to the eye - see main text). 
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FIG. 13: Contour plot of Imx - notice the "spin-pseudo- 
gap" at (g) = = Qaf = {iT/a,n/a) for energies below 
m = 20meV and the dispersion above m. 
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FIG. 12: The energy of the peak in the superconducting state 
- also called the resonance energy Eresonance as a funatjon of 
Tc in underdoped YBCO. Data taken from Dai et aB3. By 
extrapolating the linear relationship down to Tc = (dashed 
hne added as a guide) we extract an estimate for the mass 
scale related to topological gauge fluctuations 



the phase where the gauge field is in the massive phase 
due to topological fluctuations. From the data by Dai 
et a£j see Fig.|l^ we estimate the mass scale {i.e. the 
energy-gap scale) due to instantons to be of the order of 
10 — 15mey which ties in nicely with the spin- pseudo-gap 
scale ~ lOmeV estimated by Sternlieb. 

The shift of this peak to higher energies on cooling 
below Tc can be qualitatively understood as the contri- 
bution to the mass scale by condensed bosons via the 
Anderson-Higgs mechanism. This also explains the lin- 
ear relationship between Eresonance (where Eresonance de- 
notes the energy of the peak intensity in the supercon- 
ducting state) and Tc (both being proportional to a; - the 
hole doping concentration) as observed by Daiei al in 
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FIG. 14: Spin excitation spectrum for odd (open symbols) 
and even (closed circles) excitations at 5K. The open square 
indicates the energy of the maximum of the odd susceptibil- 
ity. The dotted lines correspond to the spin-wave dispersion 
relation in the insulating antiferromagnetic state with J|| 
120meV. This Fig is taken from Bourges et at 



Fig.||l2[ for underdoped YBCO. Furthermore we can ac- 
count for the very different doping dependence of the nor- 
mal state maximum (determined by the instanton scale 
with weak doping dependence) versus the doping depen- 
dence of Eresonance (determined by the condensed bosopa 
oc x) as was stressed in a recent article by Ph. BourgesEj 
From this understanding of the ASL, we can con- 
clude immediately that, Tm the temperature of the spin- 
pseudo-gap formation (and hence Tti) should be rather 
insensitive to changes in doping. Contrast this with the 
doping dependence of Tpg the formation of the pseudo- 
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gap in the single particle spectrum which as a; ^ will 
get as large as the spin-wave band width, since in that 
limit the ASL will be described by the tt-Aux phase. 

Another interesting point to discuss about the INS 
data is their momentum dependence. In Fig. [p^ we show 
the data by Bourges et al taken on YBa2Cu30Q,5{Tc — 
52) which shows the spin excitation spectrum for odd 
(acoustic) and even (optical) excitations at high ener- 
gies. The dotted lines correspond to the spin-wave dis- 
persion relation in the insulating antiferromagnetic state. 
The dispersive behavior compares nicely with Fig.]!^] 
the correlations in the ASL liquid above the mass scale 
m ~ 20meV 

So far we have argued for the strange phenomenology 
of the pseudo-gap phase to be tied to the physics embod- 
ied in the ASL. It is however clear that the ASL physics 
has to give way to the superconducting state at low tem- 
peratures. An important question is how this transition 
happens. Let us then briefly mention three plausible sce- 
narios for how the ASL goes over to the superconducting 
state. In other words, we would like to understand how 
the C/(l) gauge field gains a mass term which leads to 
the destruction of the ASL state. 

In the first picture the mass of the U{1) gauge field 
is generated via confinement due to instantons. After 
the opening of the energy gap, there is no residual un- 
broken gauge structure left at low energies. Thus the 
confinement can also be referred to as C/(l) gauge struc- 
ture breaking down to Zi. {Zi gauge structure means 
no gauge structure.) In the confinement phase {i.e. after 
the U{1) gauge field is gapped), the spinous and holons 
recombine into electrons which appear as the relevant de- 
grees of freedom at low energies. Thus we also call the 
first picture spin-charge recombination. This in particu- 
lar would suggest the existence of well defined quasiparti- 
cles at low energies (as probed by ARPES) even above the 
superconducting transition temperature Tc. (this should 
also be accompanied by resistance in d.c. conductiv- 
ity). Furthermore the biading of spinous and holons will 
gives rise to Fermi arcsJlj 

In the second picture, the U{1) gauge structure also 
breaks down to the Zi gauge structure, but now via holon 
condensation. In this case, Tm = Tc and the ASL directly 
transforms into the c?-wave superconducting state. There 
are no well defined quasi particles above Tc. 

A third scenario might be the breaking of the U{1) 
gauge structure down to a Z2 gauge structure via the con- 
densation of spipon bi-linear terms which do not break 
any symmetries. cil (a Z2 gauge structure was also found 
in a slave-ferniion description of spin systems via the con- 
densation of bosoH. bi-linear terms which break 90° ro- 
tation symmetry.E3 The Z2 gauge structure can also be 
obtained viau-the condensation of bound states of dou- 
ble vorticesJl3) The breaking from U{1) to Z2 also re- 
sults in a mass for the U{1) gauge field and leads to a 
new Z2 spin liquid, jljlit; transition to this new topo- 
logical/quantum orderclui implies the appearance of true 
spin charge separation since the Z2 gauge interaction is 




pseudo-gap 



SCspin 



e 90° e 90° 

pseudo-gap metal d-wave superconductor 



FIG. 15: Minimal quasiparticle energy Emin along the line 
in the 9 direction. The left panel depicts the pseudo-gap 
metallic state and shows a finite Fermi arc. The right panel is 
for the superconducting state. In the superconducting state, 
the quasiparticle may have four bands (instead crf-jthe usual 
two for the BCS superconductor). For details, seeLJ. 
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FIG. 16: Temperature evolution of x(cii<^) at the incom- 
mensurate peak position for energies i-gf 3 meV and 4.5 
meV. Data is taken from Yamada et ala whose sample was 
Lai,s5SrQ,i5CuOi{Tc — ZLiK). Notice the spectrum devel- 
ops a gap for the 3 meV case on cooling below Tc 

only short ranged.00Ei (Senthil experiment) 

To determine which scenario actually applies to real 
high Tc samples, we need to rely on experiments. In the 
following we would like to argue that experiments suggest 
the first scenario as most plausible in the high Tc samples. 
The appearance of Fermi arcs in the first scenario implies 
a small energy scale ^sCspin for the spin excitations as- 
sociated with the superconducting state. (See Fig.[|5j) 
Such a small energy scale was observed in the experi- 
ment of Ref. ^ (see Fig.|l^), where it was found that 
the spin susceptibility x(a;) decreases below Tc only when 
Lo < 4meV. This indicates that Ascspin is as small as a 
few meV. Direct evidence for Fermi arcs comes from-their 
observation in a recent photoemission experiment .c3 The 
expected well defined electron-like quasiparticles near the 
Fermi arcs were also observed. 



VII. CONCLUSION AND OPEN QUESTIONS 

In this paper we have given an account of the spin 
correlations in the ASL phase through the calculation of 
0{1/N) contributions depicted in Fig.||l[]. It was shown 
how the gauge fiuctuations strongly enhance the stag- 
gered spin correlations leaving the uniform spin correla- 
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tion unaffected. This result is very natural within this 
gauge fluctuation picture where the uniform spin corre- 
lation is protected by current conservation and cannot 
have any anomalous dimension. From this perspective it 
is also easy to account for the qualitative different be- 
haviors of planar Cu and O spin lattice relaxation rates 
seen in NMR experiments where C„ probes the enhanced 
staggered correlations of the ASL down to a temperature 
scale Tm ~ Tti, where the Cu 1/(2^1 T) peaks, which we 
associate with the appearance of an energy gap in the 
gauge spectrum. What is really remarkable in the ASL 
picture is that enhanced staggered spin correlations are 
obtained while at the same time having a small single 
particle density of states within the pseudo-gap phase 
without any fine tuning. As the doping x decreases, the 
pseudo-gap increases and, surprisingly, the staggered spin 
correlations also increase. This strange experimental be- 
havior is explained naturally by the ASL picture. 

It should be noted here that our ASL shows qualita- 
tive similarities with the quantum critical point scenario 
of Refs. |,|^. However as was stressed in the paper J13 the 
appearance of the pseudo-gap which they associate with 
the suppression of the spectral weight of spin waves char- 
acteristic of the quantum disordered (QD) regime should 
affect the low frequency dynamics for both q = and 
q = (7r,7r). Within this framework it then seems hard to 
account for the qualitative different behaviors seen at the 
Cu and O sites in NMR experiments. As we have stressed 
many times within the ASL this difference is protected by 
the U{1) gauge structure. Also, in our ASL approach, we 
do not assume, in contrast to the quantum-critical-point 
approach, any nearby symmetry breaking phase and we 
do not require any strongly fluctuating order parameters 



to give us critical behavior. The ASL can by itself appear 
as a stable quantum phase, or as a phase transition point 
between two states with new kinds of order - quantum 
order. The two states-.can have the same symmetry and 
no order parameters.cl Hence the ASL associated with 
the transition between quantum orders can show scaling 
properties divorced from any critical point studied so far. 

From this perspective the ASL is just one of a whole 
slew of possible quantum orders characterized by a mass- 
less U{1) gauge field coupled to massless Dirac fermions. 
It is also clear that the destruction of the ASL via the 
transition into a phase with a massive U{\) gauge field 
(phenomenologically described above via the introduc- 
tion of the mass-scale m) demands a much more careful 
analysis and highlights the main difficulty with the slave 
boson approach to the t J-model. Within this scheme it 
currently seems to be impossible or at least very challeng- 
ing to describe this energy/temperature regime around 
m/T,n theoretically. The reason for that is related to 
the need of introducing a scale which is not tied to bo- 
son condensation and hence divorced from the mean-field 
energies and the corresponding degrees of freedom. 

Through our examination of the experimental data, 
we find that the ASL, plus spin-charge recombination at 
lower energies, provides a consistent and natural (with 
no fine tuning) description of underdoped cuprate super- 
conductors. This is the main result and the bottom line 
of this paper. 
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to thank A. Seidel, for his mathematical insights. This 
work is supported by NSF Grant No. DMR-01-3156 and 
by NSF-MRSEC Grant No. DMR~98-08941. 
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APPENDIX A: SPIN OPERATORS IN THE CONTINUUM LIMIT 



In this appendix we derive the expression for the spin operators near q = (0,0), (7r,7r) and (7r,0) in terms of the 
spinon field ^ defined near the four nodes in momentum space. 
Let us look at the staggered correlation in detail. 



S(Q + q) = I E /^(P - q)^/(P + Q) 



(Al) 



where Q — (n/a.Tr/a) and we have suppressed the frequency index. After introducing 



/e^^(/(q) + /(q + Q)) 
/o^i=(/(q)-/(q + Q)) 



the expression for the staggered spin operator can be rewritten 

5(Q + q) = J E IfeiP q)> /1(P - q)] ^(^^3 - ^'J2) (j^j^J) 



(A2) 



(A3) 
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Splitting the sum J^p = Eq+k + Eq.^, where Q = (7r/2, w/2) Q = (-7r/2, -7r/2) yields 



^(Q + q) = i5^[/l(k-q+Q),/t(k-q + Q)]f(a3-ia2)(^^* 



7e(k + Q)) 
,(k + Q) 



+ 



(A4) 



Noting Q = Q + Q and using /e(p + Q) = /e(p) /o(p + Q) = -/o(p) we arrive at 

5^(Q + q) = ^E[/enk-q + Q)Jonk-q + Q)]r>3(-^^^^^^^^^^ 

= iV'^i(k-q)lf*i(k) (A5) 
k 

where ^i(k — q) = ^/j(k — q + Q), — /J(k — q + Q)^ and k is measured away from the node of type 

1 = Q = ±(7r/2, 7r/2). The remaining contribution arises from the other two nodes in the spectrum at ±(7r/2, —n/2). 
Following the same calculation for the uniform correlation function 



^(q) = jE/^(p-q)^/(p) 



(A6) 



we find 



^(q) = jE*i('^-^)'^3r*i(k) 



(A7) 



with ^'i(k — q) = ^/j(k — q + Q), — /](k — q + Q)^ and k is measured away from the node of type 

1 = Q = ±(7r/2,7r/2). As before the remaining contribution arises from the other two nodes. 

The correlation near q = (tt, 0) is obtained in a similar fashion the only difference here is that for this momentum 

transfer the two types of nodes get mixed and hence when we split up the sum over momentum space we need to 
consider all four nodes simultaneously. Under this proviso the calculation goes through as above and we end up with 



^(Qx + q) = ^E*(k-q)(^^^ ''o)^*('^) 



(A8) 



where 

*(k - q) = [/t(k - q + Qi), -/t(k - q + Qi), -/t(k - q + Q2), /t(k - q + Q2)] 
with Qi = (7r/2,7r/2) and Q2 = (7r/2,— 7r/2) 

APPENDIX B: TWO-LOOP CALCULATION OF SPIN CORRELATION 

We employ a four dimensional representation of the Dirac algebra {7^, 7y} = 25^,^ (/i, v = 0,1, 2) in the main body 
of the paper ie Trl = 4. To perform the Tr over the spinor indices we need the following identity 

which can be simply derived from the Dirac algebra by commuting 7-matrices through and using the cyclic property 
of the trace Tr. Using the above identity we can simplify 



[1(^)] = 



{2'kY J {2'kY 



Tr 



N 



{p + kYk^{k + q)'^k'^\q\^ 



(Bl) 
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to 



64 

[1(^)1 = T7 



d'^k r d'^q [k'^{p + k)-q{k + q)-q-2{p + k)-k{k + q)-qk-q\ 
N J {2t:Y J (2^ 

and for diagram |0(C)] we obtain 

d'^k r d'^q 



(27r)'* J (27r)<^ 



Tr 



(p+/c)2(A: + ^2fc4^ 

8 {p + k),Yl''{P + k + qW{k + q)0-/^rka7''{S^uf - q>,qu) 



64 



d'^k 



d-^q 



N {p + k)^{p + k + q)^{k + q)^k^\q\^ 

F{q,k,p) 



N J {2t:Y J {2t:Y {p + kY{p + k + q)^{k + q)^P\qf 



(B2) 

(B3) 
(B4) 



where 



F{q,k,p) = 



P (^fP -^k.qf + 2k-q^ + k-pq- p- ^ff 

+ q)^{2^k ■ p + fk ■ q- 2k ■ pq- p- ^<ff) 
{k + p)^{^p-q + ^k-p-\^f) 



1^-2 ^1 1-2-4 

-qp -pqq-p+-qp 



-4 

q 



In both expressions (B2,B3) the integration over k is convergent in c? = 3 and can be performed noting that 

d^k 1 1 



{2^YP{k+p)-^{k + q)^ 8|pl|<7l|<7-:^ 



d'^k 



1 



1 



(27r)3 k'^{k+p)'^{k + q)'^{k+p + q)'^ ' P 

Thus after the k integration we arrive at 



1 



1 



d'^q 



N J {2ttY 



p-q 



2p- q + p 

WW+^ \p\f\p + g] W 



\p-q\ W+p\ 



\p\ 



d'^q 



N J {2ttY 



\p\ f + 2p- q 



(ff + 2p* 



Adding the contributions 



2 X [1(A)] + 1(C)] 



d'^q \2 A\p\ p-q q^f + 2if 



N J {2ttY 



q^ q^\p + q\ \p\q^\p + q\ P ■ q(p\p\\p + q\ 



(B5) 
(B6) 

(B7) 
(B8) 

(B9) 



As noted previously because of the p ■ q term in the denominator of the last term above we wont use dimensional 
regularization. Therefore we set d — 3 and introduce an upper cut-off A to regularize the above integrals. Also note 
that we will neglect the first term in (B9) which is linearly divergent whose appearance is tight to the regularization 
via a momentum cut-off ( such divergences do not appear in Lorentz invariant regu larization schemes). 

Under this proviso we obtain the following result for the last three terms in (B9) 



4 
N 



m 4bi A2 2\p\ \p\ A2 \p\ ^^,^m, 



where the last term ( B9 ) is given by 



d^q <fp^ + 2p^ 
{2TrYp-qf\p\ \p + q\ 



87r2 



\P\ 



2 + y V^Ty-^ 
dy — In — 

y V 1 + y V^ + y + Vy 



A^ 



, 2 + y yr+^-i 

dy In — 

yVTT^ yrr^+i 



= -TfV Ci+C2 + 2ln— 



(BlO) 

(Bll) 
(B12) 
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and 



2 + v Jl + y - .fy 
Ci - / dy^=J=lvA^==L — ^ = 7.748128723 

C2 = 4?n(\/2 -1) + 4 + /n(3-2\/2)2-2;n(3- 2^2) + 2/n(3-2V2)\/2 = 2.121475665 
Finally extracting the logarithmically divergent terms we arrive at the result stated in the body of the text. 

APPENDIX C: ANALYTIC CONTINUATION OF IjN CORRECTION 

Having obtained the 1 jN correction we would now like to reexponentiate our result in the form 



This however immediately confronts us with the problem that after analytically continuing (CI) the susceptibility 
has the wrong sign for 2z^ > 1. It is however hard to understand why the spin operator cannot have an anomalous 
dimension bigger than ^ . In order to analyze this issue it is helpful to look at the corresponding Euclidean real space 
correlations. Let's take the case of a general staggered spin correlation of the form 

(5+(f)57(0)) = ^(-l)-^ (C2) 
where the spin operator has anomalous dimension v. We now Fourier transform this to obtain 

ik^'-''' r zi^^Mn) (C3) 







For < \Re{2i' — 2)| < 1 this can be evaluated^ to give 



—r(2iy-2)sin((iy-l)Tr)\q\^-^'' (C4) 
47r 



This leads after analytic continuation to the result equation! 16). From here it is easy to see how the F-function 
and the sin which are missed in the naive exponentiation of the Euclidean momentum space result conspire to give 
the correct sign for the spin correlation no matter what the size of the anomalous dimension. In particular we can 
continue the above for the case v ^ to give the correct mean-field result. 



APPENDIX D: TWO-LOOP CALCULATION IN THE MASSIVE PHASE OF THE U{1) GAUGE FIELD 



This appendix gives some details of the calculation leading to the 0{1/N) corrections to the staggered spin corre- 
lation in the case when the gauge field is in the massive phase with propagator 

8 / q^q„. 



We can use the results derived in Appendix A for the integrals over k ( P39D since they are unaffected by the change 
in the gauge propagator. 



2[1(A)] + W 



N J {2ttY 



p-q 



(fp^ + 2p* 



(Dl) 



\q\\/q^ + m^\p + q\ \p\\q\\/cp + m'^\p + q\ P ■ q\q\\/(p + ni'^\p\\p + q\ _ 
where we have neglected the linearly divergent term as discussed above. After integration over q in d — 3 we arrive at 



4 

N 



-In 



2A 



\/ p^ + 



+ 



2A 



\p\ + v/p^ + m2 



67r2 



(D2) 



2 + 2/ 



8^ I 



2 + 2/ 



VT+y-1 
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FIG. 17: Comparison between /I 
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FIG. 18: Comparison between the analytically continued forms of /I, /2 and the numeric evaluation of (D3) after analytic 
continuation (m = 30). 



It is not hard to check that in the Umit ■'^ ^ the terms proportional to m cancel. In extracting the logarithmic 
term we have to take a closer look at the last term vifhich can't be evaluated in a closed form. In Fig.(|l^) we compare 
two different functional forms with the numeric evaluation of 



\P\ 



dy 



Note that there is virtually no difference between fl 



2 + ?y yr+^-i 

= — In , 



and /2 



-A\p\ln- 



2A 



(D3) 



However this 



comparison has to be taken with a grain of salt as we are still in Euclidean space but are interested in the analytically 
continued forms. In Fig.(jl^) we plot the analytically continued {ipo — > w + ie) forms of the above functions. From 
this comparison it is now obvious that fl which has developed a pole at m is unsuitable a s an approximation and we 
take /2 which fits the numerically integrated form quite well. Collecting the Zn-terms in (D2) we obtain 



2 X [l{A)]r 



2A 



Analytically continuing and combining with the mean-field result (Wm gives 



(D4) 
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APPENDIX E: LOW ENERGY EFFECTIVE THEORY FOR THE SUPERCONDUCTING STATE 



Starting from the mean-field Hamiltonian for the spinons 



in the d-wave paired state 



al + 

we obtain after Fourier transformation the foUowing Lagrangian 

-icj + e(q) + a% 



L = 5][/|(q,c^),/i(-q,- 



-?7(q) \ ( /T(q'^) 



e(q) = ~2Jx{cos{qxa) + cos{qya)) ri{q) ^ -2JA{cos{qxa) ~ cos{qya)) 
which results with 

5+(q,w) ^ i ^/|(k,Wfc)/i(k + q,t^fe+t^q) 



S (q,w) = ^ ^/|(k,Wfe)/|(k + q,Wfe+W5) 



(El) 



(E2) 

(E3) 
(E4) 



(E5) 



(5+(q,c.)5-(-q,-c^)) = -5] 



1 -e^ (-ztJfc - £(k) - gg) (zcjfc - iujq + e(k - q) + ag) - ?7(k)7y(k - q) 



: + (,(k) + alf + riikr] [k - + (e(k - q) + a^f + r;(k - q)^ 



(E6) 



next we split = X^Qj+p where Qi with i = 1...4 corresponds to {n/2, 7r/2), (—7r/2, — 7r/2), (7r/2, — 7r/2), (— 7r/2, 7r/2) 

in this order. To extract the correlations near Qaf ~ (tt, tt) we write q — Qaf + k 

Concentrating on one of the terms in the sum - e.g. about Qi — (7r/2, n/2) we can expand 



e(Qi + p) ~ Vfpi 77(Qi + p) -V2P2 Vf = 2\f2xJ V2 = 2V2SJ pi = ^J'^ p2 



V2 



-Px + Py 



to obtain 



(5+(Qaf + k, 0^)3- (-Qaf - k, -u:)) 

1 ^ (-^^p - Vfpi - al) (iujp - iuk - Vfjpi - h) + af,) + V2P2V2{P2 - ^2) 

^ Qi,i..;p + (w/Pi + aj^)^ + vIpI (ujp - ujkf + {-Vf{pi - ki) + alf + vl{p2 - fc2) 



Qi,Wp,i=2...4 



Note ki = fc2 = Now we approximate the sums by integrals and define 



Pl = VfPi + P2 = V2P2 
ki = Vfki + 2aQ ^2 = ^2^2 



which allows us to rewrite the first term in (E8) 



(f3p p- (p-k^ + ipi 

f {v - fc) 



LOi. — iuink\ 



4 j {2n)^v-fV2 



(E7) 



(E8) 



(E9) 
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Using the Feynman trick or otherwise it is not hard to convince oneself that the last two terms cancel each other and 
the final result is 

Combining all four contributions from the four nodes in this way and translating the fc's back to the fc's which measure 
momentum along the ab axis away from Qaf we obtain after analytic continuation 

Im(5+(c^, Qaf + k)S-{-uj, -Qaf - k)) 
1 



&AVfV2 



I 






e [lo^ ~ {vfki - 


- 








vlk^^ 


9 [uj'^ — [vfk2 — 




vlkl] 



_ K + ky_ k, + ky 2V2aJx V2 = 2V2aJA 

V2 V2 



(Ell) 



the result stated in the main body of the paper. 
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